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1. Introduction 

Supersymmetry is an important guiding principle for extensions of the standard model. The 
symmetry relates fermionic and bosonic particles. Its algebra is a non-trivial extension of the 
Poincare algebra by fermionic generators. The anti-commutator of these generators is related to the 
generators of the translations. A fundamental prediction of unbroken supersymmetry is a pairing of 
the energy states in the spectrum. Consequently for each fermionic particle there exists a bosonic 
particle with the same mass. 

A supersymmetric extension of Yang-Mills theory (super Yang-Mills theory) corresponds to 
the gauge sector of supersymmetric extensions of the standard model. Fermionic fields (gluinos) A 
are added as fermionic partners of the gauge bosons. The matching of the degrees of freedom in the 
bosonic and fermionic sector requires them to be Majorana fermions in the adjoint representation 
of the gauge group. The continuum Lagrangian has the following form 



where F^y is the field strength tensor and the gauge group is SU {Nc). If the mass term nig is zero 
supersynmietry is predicted to be unbroken, even in the quantized theory [1]. A second symmetry 
of this theory is the (chiral) U{\)r symmetry, X e~^^^^X. It is broken by the anomaly. The 
remaining Z2Af^ is spontaneously broken to Z2. 

The theory is assumed to be confining, therefore the low-energy degrees of freedom are formed 
from gauge invariant compound operators of the fields. The symmetries of the theory constrain the 
form of the low-energy effective action. The first prediction for a low-energy effective action was 
derived from one supersymmetry multiplet that contains the meson operators a — fo (XX), a — 7]^ 
(A75A), and the fermionic gluino-glue (cTuv^/iv^^) [2]. Later on, as an extension of the low energy 
effective theory, a multiplet of glueball operators (0++ and 0~+) and the fermionic gluino-glue was 
added [3]. 

2. Super Yang-Mills theory on the lattice 

Our approach is based on the work of Veneziano and Curci [4]. The lattice action consists 
of the well-known plaquette action together with the Wilson-Dirac operator in the fermionic 
sector [5], 



(1.1) 




(2.1) 



The Wilson-Dirac operator is defined in the following way 
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(2.2) 



where K = 



2K+4) 



and the lead to a gauge invariant shift of the fermion fields in ji direction 



r^A(x) = y^A(x+A)- 



(2.3) 



2 



Simulations of super symmetric Yang-Mills theory 



Georg Bergner 



The links in the adjoint representation, V^, are generated from the fundamental links using {y^)ah = 

2Tr[c/^r^c/^r^]. 

Supersymmetry and the U{\)r symmetry are broken in this discretization. Supersymmetry 
is known to be broken in any local interacting lattice theory [6]. The U{\)r symmetry is, like 
chiral symmetry in QCD, broken by the Wilson mass. The breaking is not introduced in terms of 
a Ginsparg-Wilson relation, that ensures a modified symmetry on the lattice. Hence a fine tuning 
of the parameters is needed to recover the symmetry in the continuum limit. Fortunately it turns 
out that the number of terms, which must be considered in the fine tuning, is small - in case of the 
lattice action 2.1 actually one. 

The Ward identities of supersymmetry and the C/(l)/? symmetry on the lattice are 

(V^/,^(x) ff{y)) = m,{Ds{x) ^(j)) + {Xs(x) (2.4) 
(V^/^(x) ff{y)) = mg{DA{x) ff{y)) + {Xa{x) i^{y))+ o< {FF ff{y)). (2.5) 

The expectation values on the left hand side contain derivatives of the super-current Js and the 
current of the C/(l)/? symmetry Ja with an arbitrary operator composed from the fields of the theory. 
A conserved symmetry requires the right hand side of these equations to be zero. The terms Ds 
and Da also appear in the continuum and represent the symmetry breaking induced by the non-zero 
gluino mass nig. The Xs and Xa{x) are induced by the discretization and vanish in the classical 
continuum limit (tree level). The last term in Eq. 2.5 stands for an additional contribution of the 
anomaly and is present also in the continuum. 

In a quantized theory one has to include operators that can mix with Xs and Xa in the continuum 
limit. The analysis of these operators for chiral symmetry and supersynmietry can be found in [7] 
and [4]. In our case it leads to the following expressions 

{V^ZAJ^{x)ff{y)) = {mg-mg){DA{x)ff{y))+oc{FFff{y)) + 0{a) (2.6) 
{V^{ZsJ^{x)-^ZsJY{x))^{y)) = {mg-fng){Ds{x)^{y)) + 0{a). (2.7) 

This means that up to effects that vanish in the continuum limit, the symmetry breaking of the dis- 
cretization can be absorbed into a redefinition of the currents and a shift of the bare mass parameter 
rUg. A tuning of rUg is hence enough to restore the chiral symmetry and supersymmetry in this 
theory. Up to 0{a) effects the tuning for both symmetries is the same. 

The above considerations show that the chiral and the supersymmetric Ward identities can in 
principle both be used to tune the bare mass parameter. However, the signal of the supersynmietric 
Ward-identities does not allow for a precise tuning. The chiral Ward-identities, on the other hand, 
contain an additional term from the anomaly. To avoid these limitations, we consider the con- 
nected part (OZI approximation) of the a — rj', the adjoint pion (a — k). The mass of this particle 
vanishes in the chiral limit. However, \ht a — K can be defined only in a partially quenched setup. 
Therefore, the consistency with the supersynmietric Ward identities must be verified in a separate 
measurement. 

3. Some details about the simulations 

The simulations of this theory include several specific challenges. The integration of the Ma- 
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Figure 1: The left figure shows how our method focuses the calculation on the real eigenvalues. The 
red boxed points are the part of the complete spectrum of computed in the iteration. The 20 lowest 
eigenvalues in the right figure are obtained from 2007 configurations at /c = 0.1495, j8 = 1.75 on a 32^ x 64 
lattice. Only 16 negative Pfaffians have been obtained. 

jorana fermions leads to the Pfaffian Pf(CDw), where C is the charge conjugation matrix, instead 
of the determinant. The modulus of the Pfaffian is the square root of the determinant. It can be in- 
cluded in the update algorithm using a polynomial approximation (PHMC algorithm [8]). The sign 
of the Pfaffian, however, can only be included in a reweighting of the measurements. A detailed 
analysis shows that the Pfaffian has the same sign as the product of the doubly degenerate real 
eigenvalues [9]. If the number of negative real degenerate eigenvalues is odd, a negative Pfaffian 
sign must be included in the reweighting. Our efficient way to calculate the real negative eigenval- 
ues using the Arnoldi algorithm together with polynomial focusing and acceleration, illustrated in 
Fig. 1, can be found in [10]. 

The considered observables of the theory are the correlation functions of glueballs, mesonic 
operators, and the fermionic gluino-glue. A measurement of the mass of the glueball is possi- 
ble using variational smearing methods (see Fig. 2). 
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Figure 2: The effective mass of the glueball obtained with several different methods. The variational smear- 
ing methods allow the determination of a plateau of the effective mass. (This result has been obtained on a 
32^ X 64 lattice, j3 = 1.75, k = 0.1492, 2848 configurations.) 

All of the meson operators in this model contain disconnected contributions, that become 
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dominant in the chiral limit. There are several techniques to obtain these contributions like the 
stochastic estimator technique (SET), or the improved volume source technique (IVST) [11]. In 
our latest studies we have applied the SET (see [12] for a review of this method), since it can 
be combined with the exact contribution of the lowest eigenmodes. In this method the inverse is 
approximated using random vectors (G Z4). We have improved the SET using the truncated solver 
method and the exact contribution of the lowest eigenmodes. 

An approximation of the inverse is obtained from the lowest eigenmodes of the Hermitian 
matrix 75 D^;. When the noise is projected orthogonal to the eigenspace, this approximation can 
be combined with the SET. In our case we have used the eigenvalues from the even-odd precon- 
ditioned Hermitian matrix. They have been obtained with the Arnoldi algorithm using Chebyshev 
polynomials to accelerate the algorithm. The spectral decomposition using the eigenvalues of the 
even-odd preconditioned matrix combined with the SET and the truncated solver method leads to 
a considerable improvement for the determination of the <2 — T]' at small gluino masses. Note that 
the same eigenvalues are used to calculate reweighting factors to improve the polynomial approxi- 
mation in the PHMC algorithm. An illustration of these methods can be found in Fig. 3. 
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Figure 3: This picture shows different strategies to obtain the disconnected contributions of the a — if]'. 
Different number of stochastic estimators in the SET are compared with a truncated eigenmode expansion 
using 10 and 100 eigenvalues. The last two strategies shown in this picture are a combination of the SET 
and the truncated eigenmode expansion. In the last step a truncated solver method (using 100 low precision 
and 10 precise estimators) has been applied. The projection of the noise to the space orthogonal to the 
lowest eigenmodes leads to a faster conjugate gradient inversion. Therefore, the time needed for the last 
three strategies is comparable at small gluino masses. (The disconnected correlator in this figure has been 
obtained on a 32^ x 64 lattice, j3 = 1.75, k = 0.1494, on 300 independent configurations.) 

The fermionic particle considered in the mass spectrum is realized by the gluino-glue operator 
<7^^Tr[F^v>^], where F^y is represented by clover plaquettes. A combination of APE smearing on 
gauge fields and Jacobi smearing on A has been used in this case. This combined smearing leads 
to a good signal compared to glueballs (see Fig. 4). 

Our simulations are done with the PHMC algorithm with determinant breakup. A tree level 
Symanzik improvement of the gauge action and stout smearing for the fermion links has been added 
to the standard setup presented in the previous section. It is expected that these improvements 
reduce the lattice artifacts. 
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Figure 4: The effective mass of the gluino-glue. The correlator of the gluino-glue has been obtained with 
APE and Jacobi smearing. The gray Hne indicates the value obtained in a fit. (This result has been obtained 
on a 32^ x 64 lattice, j3 = 1.75, k = 0.1492, 2848 configurations.) 

4. Results 

In our first simulations we have used the lattice sizes 16^ X 32, 24^ X 48, and 32^ x 64. We 
have chosen j8 = 1.6 and measured the Sommer scale (ro). If one sets ro = 0.5fm, like in QCD, the 
lattice spacing is below 0.088fm and the lattice size is around 1.5 — 2.3fm. The smallest adjoint 
pion mass has been at about 440MeV. The tuning of the bare gluino mass towards the chiral limit 
was found to be consistent with the supersymmetric Ward identities. Details can be found in [13]. 

The mass spectrum obtained in these simulations revealed an unexpected mass gap between the 
bosonic and fermionic particles. Such a mass gap is in contradiction to supersymmetry. Therefore 
a more careful analysis is essential. 

As a next step of our scrutiny of the mass spectrum of supersymmetric Yang-Mills theory we 
have investigated the influence of the finite lattice spacing. To decrease the lattice spacing we have 
increased /3 from 1.6 to 1.75. The simulations have been performed on a 32^ x 64 lattice. The 
lattice spacing in these simulations is around 0.057fm with a lattice size of approximately 1.8fm. 

The preliminary results we have obtained with these new parameters indicate a major influence 
of the lattice artifacts (see Fig. 5). The most evident effect is the reduced mass of the gluino-glue. 
Hence a much smaller mass gap is found in the spectrum. A naive linear extrapolation based on 
the difference between the gluino-glue and the a — T]' towards the continuum limit leads even to 
a negative mass gap at zero lattice spacing. For a precise investigation of the continuum limit we 
need, however, to obtain more data and investigate further improvements of the action. 

5. Conclusions 

Theoretical considerations indicate that supersymmetric Yang-Mills theory with Wilson fer- 
mions on the lattice can be tuned to a supersymmetric chiral continuum limit. The simulations of 
this approach are challenging and advanced techniques must be used to get a reasonable signal of 
the observables. Besides the technical difficulties the theoretical considerations do not quantify 
the magnitude of the lattice artifacts when the continuum limit is approached. Our new simula- 
tions indicate a major influence of these artifacts on the spectrum of the theory. In particular, the 
mass gap between bosonic and fermionic particles at a larger lattice spacing is inconsistent with 
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Figure 5: The mass spectrum at a smaller lattice spacing. In the left figure the masses of the <2 — r]' and the 
gluino-glue are extrapolated to the continuum limit. The values obtained at a smaller value of j8 are included 
for comparison. The arrows indicate the mass gap in the chiral limit. The right figure shows the masses of 
all particles obtained at j3 = 1.75. As expected the data suggest a mixing of the a — and the glueball. 

supersymmetry. This gap might, however, vanish in the continuum limit. Further simulations and 
a further refinement of the approach are needed to prove these findings. 
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